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Abstract. In this paper a linear and a non-linear fractional differential equations are considered 
to solve the two differential equations analytically by Power Series method (PSM), Adomian 
Decomposition Method (ADM) and Homotopy Analysis Method (HAM). Solutions obtained by above 
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1. Introduction. Fractional differential equations are generalization of classical 
differential equations of integer order and have been used successfully to modify many 
existing models of physical process for improving results. Fractional calculus is a 
field of mathematical study that grows out of the traditional definition of calculus, 
both integral and differential operators in much the same way as fractional exponent 
is an outgrowth of the exponents with integer value. Apart from diverse areas of 
mathematics, fractional differential equations arise in rheology (Trmohoy, 2003), visco- 
elasticity (Soczkiewicz, 2002), Ultrasonic Wave propagation in human cancellous bone 
(Sebaa, Fellah, Lauriks and Depollier, 2006), Modelling of speech signals (Assaleh 
and Ahmad, 2007), Modelling of the Cardiac Tissue (Magin, 2008), Sound waves 
propagation in Rigid Porous Materials (Fellah and Depollier, 2002), Lateral and 
Longitudinal control of autonomous Vehicles (Suarez, Vinagre, Calderon, Monje and 
Chen, 2004), Fluid Mechanics (Kulish and Lage, 2002), Edge detection (Mathieu, 
Melchior, Oustaloup and Ceyral, 2003) and in many other areas. There is no doubt that 
the fractional differential equation has become an exciting new mathematical model for 
representing many diverse problems in mathematics, science and engineering. But the 
solutions of these equations involve rigorous mathematical techniques. 
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The concept of fractional calculus (fractional derivatives and fractional integral) 
is not new. In 1695 L’Hospital asked Leibnitz, one of the inventors of calculus the 
d"y 
da”? 


Leibnitz replied that “d'/?x will be equal to zv/dr : a”. After that discussion on 
this idea was untouched till the beginning of 20th Century. 


question as to the meaning of if n — 1/2 i.e., "if n is fraction"? 


The whole theory of fractional derivatives and integrals was established in the 
2nd half of the 19th Century. Caputo and Mainnardi (Caputo and Mainnardi, 1971 
and Caputo, 1966) and Caputo (1974) found good agreement with experimental 
results when using fractional derivatives for description of viscoelastic materials and 
established the connection between fractional derivatives and the theory of linear 
viscoelasticity. Adolfsson et al (Mainardi, 2010) constructed a new fractional order 
model of viscoelasticity. 


A lot of works have been done by using fractional derivatives for a better description 
of considered material properties. Several numerical methods for solving fractional 
differential equations have been evolved. À new analytical approach that can be applied 
to solve nonlinear fractional differential equations (NFDE) is to employ the Homotopy 
Analysis Method (HAM) (Golmankhaneh, Golmankhaneh and Baleanu, 2011 and Liao, 
1995, 1997, 2003). The HAM solution for Nonlinear Fractional Differential Equation 
(NFDE) were presented by D. Das (Das, Ray, Bera and Sarkar, 2015). Adomain 
Decomposition Method (Adomian, 1988, 1994) provided analytical solutions as well 
as numerical approximate solutions to both linear and non-linear differential equations 
without linearization or discretization. The ADM is also implemented to give analytical 
solutions for the class of Lane-Emden Equation (Das, Ghosh, Datta, Bera, Sahu and 
Ray, 2015). These methods are different from that of perturbation method. 


Recent seminars, symposia and proceedings have highlighted the application of 
fractional calculus in physics, continuum mechanics etc. 


So, it has become necessary to solve fractional differential equations. Several 
methods are available to solve fractional differential equations, for example PSM, ADM, 
HAM, ete. In this paper we solve a linear and a non-linear fractional differential 
equation and have shown that the solutions obtained by the above methods are found 
to be exactly same calculating up to nine terms. If we proceed further, it can be shown 
that the solutions by the above methods will be the same up to the desired terms. 


2. Formulation of the Problem, for linear fractional differential equation 


2.1 Mathematical preliminaries on Fractional Calculus. To solve the 
problems of fractional differential equations many definitions of fractional calculus 
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have been developed. The most frequently encountered definitions include Riemann- 
Liouville, Caputo, Wely. Rize fractional operator. We introduce the following 
definitions [23,24] in the present analysis. 


2.1.1 Definition. Let o € R^. The integral operator I^ defined on the usual 
Lebesgue space L(a,b) by 


rro = AS = f emn 


I f(x) = f(x) 


for x € [a,b] is called Riemann-Liouville fractional integral operator of order a (> 0). 
It has the following properties : 


(i) I? f(x) exists for any z € [a,b]. 
(ii) II? f(x) = I**? f(a) 
Gn) JT f(x) = TP TF 


T(y +1) 


G 
DO ee 


x°t? where f(x) € Lla,b], a, B > 0, y > —1. 


2.1.2 Definition. The Riemann-Liouville definition of fraction order derivative is 


RL pa = d” n—a EN 1 d” 7 n—-a—1 
DEN = q LR) = gua L 6-077 10d 


where n is an integer that satisfies a € (n — 1,n). 
2.1.3 Definition. A modified fractional differential operator ¿DF proposed by 
Caputo is given by 


cpa — n—a d" m 1 ii n—a-1 ein) 


where a € R* is the order of operation and n is an integer that satisfies a € (n — 1,n). 


It has the following two basic properties [25] : 


(i) If f € Los(a, b) or f € Cla,b] and a > 0 then 5D2 oI2? f(x) = f(x) 
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(ii) If f € C”la,b] and if a > 0 then 


n—1l p(k) 
RHA = f) EI 


2.1.4 Definition. For m being the smallest integer that exceeds a, the Caputo 
time fractional derivative operator of order o > 0, is defined as in [26] 


ulw, £) T 
Dřu(x, t) = u(x,t) _ SS? deep 
a 7 ot? 7 1 t Out E) 
m=0a—=1 y ; _ 
mcs £) —3d di m-lSacm 


2.1.5 Relation between Caputo derivative and Riemann Liouville deriva- 
tive: 


m-l p(k) 
CO _RLpoa f (0+) k-a 
oDr f(x) = o Di f(x) PEST 


; a € (m—1,m) 


Integrating by parts, we get the following formulae as given in [27]: 


O Taste f ert ateiger [Fn np f] 


j=0 


a 


b 


(ii) For n= 1, [ se )E Fas dée f f(x)®2 Dp aL )dz + Lu “g(a Inte 


a 


2.2 Linear fractional differential equation. We consider the linear fractional 
differential equation 


dy | *d^y(z) 
ind = 2.1 
Y (a) =e (2.1 
where a, 0 <a < 1 is a rational number and initial condition is y(0) = 1. 
d*y(z) — cpa : 1 1 fy? 
Here de. 7 D*'y(x) is the Caputo derivative of y(x) of order ‘a’. 
x 
d ege 
We define D = q æ nd ^D" = qa then the differential equation (2.1) takes the 
z 


form 


Dy(a) +“D*y(x) + y(x) = x (2.2) 
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We consider 0 < a < p/k < 1, where p and k are integers which are prime to each 
other. 


Then p< k > k-—p>0O. 

3. Solution of the problem by different methods 
(a) Power Series Method: 

Applying 1% on both sides of equation (2.2), we get 


I*Dy(a) + 1° *D^y(z) + I%y(x) = I"(z) (3.1) 


The integral operator I“ according to the Riemann-Liouville fractional integral operator 
of order ‘a’ is: 


Iex = Iti ate (3.2) 
a+y+l 
S v (0) 
and I*°D°y(x) = y(x)— 5 d x1, where m — 1 «€ o « m. In this paper, we have 
q=0 ` 
0O<a<l, ie, ,m=1 
I**D'y(z) = y(x) — y(0) (3.3) 
Let 
S p 
y(x) = 5» anx”/*, where 0<a<?<l1 (3.4) 
n=0 k 
ft n 
Then Dy(x) = 5 (2) De 
n=1 
Om fy [Z re ltk 09 cn E " 7 
pte) = È (Pap ES (alos 
n=1 k Ce a n=1 k E 
29 ZU n-k4p 
= Y Oe m 4 (3.5) 
n=1 k 
oo 
and I**D*y(z) = ylz) — y(0) = 3 ans” (3.6) 
n=1 
with ay = 1, from (3.4). 
Again 
Q = Eer 
ee EE 2 dE Y + (3.7) 
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Substitute the values of IY Dy(x), I^*D*y(r) and I*y(x) from (3.5), (3.6) and (3.7) 


respectively to the differential equation (3.1) we get 


xs (n ` k) uL Dale L X" aan 


n=1 


Do e a 
— Tm + p + k)/k 2+0 


From the last term of L.H.S. of (3.8), 
an =0 for n=1,2,---,k—p. 


Putting j = n + p in the third term and rewriting od we get 


|-k—p»/k ; (p+k)/k 
E EDE E: 


LT GrS/k CFOM 


Putting j = n — (k — p) in the 1st term, then 


S (+ (2k — p))/k 
zd CET Ma Ei 


= T Ka SES 2 + (p/k) 


From 1st and 2nd term, for j = 1,2,---,p—1, 


añ EE EE 
"TT UTK 


+aj=0 


(3.8) 


(3.9) 


(3.9a) 


(3.10) 


From (3.9) and (3.10), we get for =1,2,*--,k—1, aj; =Oie, k=2, j 


From (3.10), we can write 


=> L + (2k — p))/k ilk j/k 
CENT P 


+Ya j-p Ges 


= KE [GAR ^ — Rr 
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For j > p (and j 4 k + p) 


U + (2k — p))/k (G+k—p)/k 
Dis +aj;+a 0 (3.12) 
TE | TT [i8 
For = k + p, 
E [2 1 
42% T + Aktp + ak = (3.13) 
TO 
From (3.12), j = p 
[2 [1 
P i AE (3.13a) 
(k+p)/k ^  [(k+p)/k 
=> ak = ag as ay = 0 from (3.11) 
y(0) = ag = 1 (given) > a, =-—1 (3.14) 
For j = p + i (from 3.12) (j Æ k) andi = 1,2,3,--- 
(2k + i)/k (k 4- i)/k 
apa A ge A 3.15 
CT [(k 4- p 4 i)/k SS 
4+% j 
From (3.15), A244" I 0/2 di 2+ E zi 
(3+1)/2 3 + i)/2 
For i = 1,3,4,5, (i £ 2). 
[5/2 1 
SEENEN €— 
a3 [2 02 a3 5/2 
[3 [2 1 1 
BER "Es “a "2 
, [7/2 [5/2 2 
For i = 3 (from (3.15)): as 44 Ie Pes E 
KS (from (3.15)): as 3 + a4 + M3 3 a5 WE 
la ed =O S 
= 4 (from (3. : d6—— +05 +a = ag = = 
6 7/2 5 4 7/2 6 6 
l [9/2 "E 1 
For i = 5 (from (3.15): az — =0 == 
or 2 (from ( ): az [4 + ag + 05 n a7 Op 
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[5 [4 2 1 
For i = 6 (from (3.15)): ag —— + a7 + ag —= = U a —— === 
( (3.15)) dercher ge 1 
11/2 9/2 1 
For i = 7 (f 3.15)): =—+ag+a == 
or à (from (3.15)): ag 5 as + a7 17 
6 1 
For i — 8 (from (3.15)): vun + ag + D a10 IG and so 
on. 
yeas, E 
For a = 5 
"M , 53/2 , jd gre " 25 r x72 yg " 9/2 " Rie 
Ste € T . 
Upsm 5/2 2! 7/2 3! 9/2 4! 11/2 


(b) Adomian Decomposition Method. We consider the differential equation of 
the form 
L(y(z)) + N(y(a)) = f) (3.17) 
where L is the (highest order) linear differential operator and N is a (in general) a 
nonlinear operator (which includes all other operators in differential equation other 
than L). 
The function f(x) is a known function of independent variable x and y(x) is the 
unknown function which is the solution of (3.17). 


ADM defines the unknown function y(x) by an infinite series 
oo 
y(z) = $ yilz) (3.18) 
i=0 


where the components y;(r) are determined using recursive relation. The non-linear 
operator N (y(x)) can be decomposed in an infinite series of Adomian polynomials given 
by 


N(y) = 3. Ai(yo(x), 2izl,: Le) (3.19) 
i=0 
where A;'s are called so-called Adomian polynomials of yo(x), y;(1),---,yi(x) and 
defined by 
1 [ d Gei . 
Ai = = -IN ys fori =0,1,2,--- 
i! | dX = 
2-9 A=0 
1/ d 2 3 
= = (= [N (yo + Am + ENEE (3.20) 
i! | dX T 
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Now, we consider the fractional ordinary differential equation : 


dy , “a%y(x) 


dx dz? 


c Aa 


where 0 < a < 1 and ENT 
defined by the Caputo [12]. 


+y(1)=x, zU (8.21) 


= ^D? is the fractional differential operators of order ‘a’ as 


y(0) = 1 is the initial condition of equation (3.21). 
Q 


d d 
Define D = — ep? = — 
efine d. and EES 


, then equation (3.16) can be written as 
Duel + *D^y(z) + y(2) = 2 (3.22) 
In equation (3.22), we choose L — D and 
N(y(z)) S*D?y(x) + y(x) and jf(x)-z. (3.222) 


Then rewriting equation (3.22) in the form L[y(z)] = —N(y(x)) + x where L^! = I. 
Applying L-! on both sides of equation (3.22), we get 


2 


y(x) — (0) = E 


SCT L~[N(y(2))] (3.23) 


Using (3.18) and (3.19), equation (3.23) can be written as 
y(z) = Y w(z) = yola) E" E Ai(yo(z), y(r) +++ wc (3.24) 
i=0 i=0 


From equation (3.24), 


yi(z) = —L~* (Ai-1(yo(2), mie), «+, y (2), (3.25) 


1 gp Sa 
Aja = Gop | a N | My; for i = 1,2,3,- 
l pu A=0 


where 
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For i = 1, 
yi(x) = —L~*(Ao(yo)) = L^! [FD*yo(x) + yo(x)] 
= —I'~* [yo(z) — yo(0)] — Iyo(x) 
q3-9 r3’ 
a AE 
For a = 5, 
pe e 
nl) =-2- 7777 a 
From (3.25), 
yla) = - (017? + Hun (2) 
q 2a Ain q ha T^ r2 
= 2 3.25 
b-a Ba mem di 5 Ke 
For a= $, 
go? g2 9 git gå 
= 2 
volt) = Tm s sme a 
y3(x) = “(159 + insta), ya(0)=0 
go 3 go 20 z372 oa q3-2a x3 x5 
= 3 2 2 
= S [6 — 2a T 4-a T [6 — a T [4 — 2a N 3! N a 
For a= 5, 
ie ge at. Be 4 0/9 ző 
= 2 3 3 
ya) E tratara a mg =| 
ya(z) UI" + Dys(z), y3(0) = 0 
q6-4o q6-3a q6-2a qa 
= 4 6 4 
17 — 4a T [7 — 3a H [7 Ze * T—a 
"m q 20 " gto " qp 4-38 " qt * zê 
[5 — 2a Ba [|5—3.q A4 6! 
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For a = 7 
ba x5/2 E 32? n 34772 " 2.7 " 44:9/? D 44:1/2 E 62? n go 
0 6 
œ 72 3 p2 4 fp nsa 5 6 


ys(zx) = (159 + Dyalz) 
| q Do q Aa Det q 20 gia 


5 10 10 5——— 
[Boba [S—4a =a ea T 8-a 


TE up c ud + LEG 
[6 — 3a [6 — 2a [6 0 6—4da 5! 


5—3a 5—2a qa wo 4a r5 d 


_1 
For a= 5, 


de? x’ Gil All Geh  10xH/2 1049 5 x!3/2 q 
= Plu a 4 + A 
ys(z) 92 3 4 ny s EY Ba n 


yola) = (1794 Dys(x) 
8—6a 8—5a 8—4a 8—2a 8-20 6—5a 
ZE T LE M T x 
= 6 15 20 15 
w= o ea Pena TES 
6—4a q6-3a q9-2o 6—Q q8-o xê T 


M dd 
10 10 5 5 
fms ea me^ Bow ^B s 


YT mus 
[7 — 4a 


|= 


For a= 5, 


(2) a » 5a? " 1052/2 " 112? " 1114112 16049 204137 
m. = 
SR 92° 4 ny 5  Js2 6 "Em 


157" DIE zë 
tm" mms s 


Hetz) = (179 + Dys(z) 


9—Ta 6a HBO 40 93a 


SE x 
7 21 35 35 
Ho—7a T RS ^ locom * Moda * 10—3a 


9-20 q 8a q Do q 4 q 30 
20 6 15 20 
i U- h-a a-pa" SE 8 — 3a 
q 20 q 0 al Bu q? 


15 6 6 
M Sh == m-a ny 
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For a = 1, 
4 9/2 5 11/2 
o x 15x x 
x 6 1 
y) = -7 1/2 5! 13/2 
( ) = 10-80 " qp10-7a ae gl0-Ta 17 q10-6a 
me Nr fiw (eve Mt 
n E 
For a= 5, 
79/2 
t) = SR 
ys (2) 1/2 
4/9 z2 "M if gt 0/9 
R = l= F 2 PaT MAA 
vapula) 52' 7 “F2 s p2 4 ny 


d 
c) Homotopy Analysis Method (HAM). Let us consider L = a = D and 
e 


T2 D=, 


The differential equation can be written as N(y(x)) = Dy(1)+"D“y(1)+y(1)-x = 0 


Then 


N (o(z, q)) = Dó(z, q) zb "DU" AL, q) 2x OLT, q) =x=0 


Deformation equations : 
L [y (x) m 


where 


dml 


(3.27) 


(3.28) 
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Expanding oz, q) in terms of q: 


iz, ol = dëi, OI + Y (55), E d 


m=1 oq” =0 m! 
= yo(z) + A. yn (z)a" = 3. yn (z)q" (3.29) 
m=1 m=0 
where 1 [and 
Sm (at) SC GF for m = 0,1, 2, 
For q = 1, p(x,1) = y(x) = Y ym(2) 
y=0 
For HAM, 
From (3.29), 
ort 
bj (2-5) 
(mI 7 I — 1)! gni = 
Then 


= 1 J” pa, ail a 1 a” pa, all 
Em = ( 1)! üqm-1 )) B (a 1)! | üqm-1 ). 


1o Vëloe, a) 
d Igo )) nen 


= Dym-1 +“D“Ym-1 + Ym-1 — (1 — Am)zx for m = 1,2,3,- (3.31) 
From (3.28): 
Llym(2) — Anym-1(z)] = h [Dym-1 + “D*Ym-1 + Ym-1 — (1 — Am) 2] (3.32) 
for m = 1,2,3,--- 
For m= 1, 


Llyi(x) — Myo(x)] = h [Dyo(x) + *D*yo(x) + yo(x) — x] (3.322) 
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Let yo(x) = 1 
Huel = hlyo(x) — 2] = h(1 — v) 
y (1) = hI [yo(a) — a] = h(a — 2*/2) 
For m = 2, we get from 6. 32), 
L[ys(a) — Aegm tel = (1+ h)u (a) +AU + D(a — 27/2) 


=(1 h? 
ERT E J-a 8 


For h= -1,a-1l, 


n) = a m: me s 
via) = (1 + h)yo(x) + A(T™® + mie) 


3—2a 3—a 3—a 3 
x £ x x 
—(1+h)[------ hi 
ee |+ E 
q 2a gto q a at 


5-20 [5-a [5-a A 


2 5/2 7/2 4 
x zx x x 
x)= 2 +2 + 
se) = ar "mee "Tee" 
4—3a 4-20 4-20 4—a Jo 
x T x T T 
=(1+AM..... n? 2 9 
ya = (ERO [se (eet tet Boda b-a B-a 
xt go 30 go 2a q 20 go q9»-a x’ 
+ 2 2 
4 [6—3a 6 — 2a [6 — 2a 6-a [6-a 5! 


(x) A + E +2 DE 234 3459/2 ` ch 
T m 
S 73 ^3 "^a 4 mnm 5 
5-4a x57 3a 5-20 5—a 5 
Ra HH x HH 
= (14 Ae h? 6 4 
AS a (eae 5 
q6-4o q6—3a 8-20 qó=0 xê 


4 6 4 
[7 — la I" — 3a [7 — 2a T-a 6! 
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For h=-—1,a=3, 


q3 mr m z5 gil/2 gô 
= 4 5 5 4 
ys(v) 3! 9/2 4 uu x 13/2 "e 
6-50 6—4o. 6-30 620 
EE wl = p 10.2 ip 
ye(z) = (1+h)| |+ — meum * ma" [7 — 2a 
Én H q 40 gTa 
5 5 10 
T 7—a' 6 [8 — 5a [8 — da [8 — 3a 
- q 20 qa q 
[8 — 2a B-a 7! 
"TN «il 
Forh=-—1,0= 5, 
r) = +5=— +9 +5 5 9 5 
we) = pa a ae toa ma ^s BA 7 
7—6a 7-50 T=40 
x x x 
—(1+h)[------ he 6 15 
yr(z) = (1+ h)[ E Bm "mmm "t accum 
q 30 gl 2a 
S 8 —3a N 8 — 2a = 
For h = —1, a = 5 
TM q 19/2 g^ 
Se A fim 5! 
3 Ta 360 q ¿60 zë Do 
ie EE hs 6 6 SCH 
IMA a eta mea” 
ys( 
"m ai 
Forh=-—1,0= 3, 
19/2 zê 
= +75 + 
3/2 z2 5/2 g3 g2 yê r2 
+ dk des (8.33) 


E 
e legy + 2 + 
Vil T 5/2 2! 7/2 3! 9/2 4! 11/2 
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4. Comparison of three methods used in solving a nonlinear fractional 
differential Equation: 
(A) Power Series Method (PSM): Let us consider a non-linear fractional 
differential equation, 
d 
i, —“DPy-y?=2x, 2>0 and y(0)=1 (4.1) 
Let 


y(x) = > apar? (4.2) 
n=0 


be the power series solution of the fractional ordinary non-linear differential equation 
(4.1). 
Applying 11? on both sides of (4.1), 


I? py M (pe ept/2)y u por 2 — 3/2 (4.3) 


Here ^D? is the Caputo derivative of order 1/2 and 117? is the R-L integral of order 
1/2, 


m-—1 k 
We know that 11/2cp1/2f(z) = f(z) — Y POF x>0andm-1<a<mis 
k=0 : 


the order of fractional derivative. 

Here a = 1/2 and m=1, 

Now, if f(x) =1, then SG DMS E =1-1=0 

Again, if f(x) is a function of x, such that f(0) = 0, then 112 °D"? f(x) = f(x) — 
f(0) = f(x). 

Substituting (4.2) in (4.3), we get, 


iy Y — asc ^ — ag 
ad 2 n/2 4- 1/2 


n=0 
Ro m n—1-i i 2 
Iph 2 2. On-1 047 2 22 = 5/2 13/2 (4.4) 


In the last term, put n = 7 +1 (n= 1,j = 0) and in the other two term replacing n by 
j (after integration) 


SS 12) Gel. "NES Am 
A ane Zon 


=0 


La 
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E» Gum, € 
22x63) 5/2 
d) ea cm as =1 
For j = 0, a 0 a, =0 
For j = 1, m ag — 01 en a2 02 1 
For j = 2, 23/2 En a —0 dE 
For j = 3, from (4.5), = a4 — a3 » E 43.1 idi N 04 == 
In general for 7 > 4, 
(2 5 d af Qj+1 — Qj Vot =0 
Forj=4, as = Es a4 + Lo > gp Dus 
For j =5, ag— um as + ACE = = 
For j = 6, BEER GË iai) = 0/7 
For j=7, ag = ake x 3 
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` General power series solution is: 


56 19 608 
=j E BaS 8/9 L. 7.5 7/2 
YPSM = tots ee +, ER + 6 ras 


nce 3 G 1280 102 LT 
12 ' 97)" * VI35/7s 28357 


(B) Adomian Decomposition Method (ADM). We consider the non-linear 
fractional differential equation : 


L[y(x)] + N[y(x)] = f(x) 


where Lis a linear operator and N is a non-linear operator and f(x) is a known function 
of x. 


ADM defines the unknown function y(x) by an infinite series : 


ge Kë yi(x) (4.6) 
1=0 


where the components y;(x) are determined by using recurrence relations. 


The non-linear part N[y(x)] can be decomposed in an infinite series of Adomian 
Polynomials given by 


= GËT (4.7) 


where A;'s are the so called Adomian Polynomials and defined by 


Ai(yo, yi: Yi) = ` IE: E » vun) | $—12,9,-- (4.8) 
A-0 


j=0 
Explicitly (4.8) can be written as, for N(y(x)) = “D'Py(a) + (y(x))? 


Ao = N(yo) = ye + ^D? yo(z) 
Ai = yN’ (yo) = 2yoy + *D'Pyi (2) (4.8a) 
Az = yaN (yo) + ty? N" (yo) = 2yoya + yt + *D'? y; ' 

( 


An = yaN (yo) + y1y2N" (yo) + Sf IN" (yo) = 2yoys + 2y1y2 + “DP ya 
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We will apply ADM for our problems : 
Dy(z) - °D'Py(x) — y (e) 22z, 220 


with initial condition 
y(0)=1 
Rewriting equation (4.9), 
L[y(z)] = Nly(x)] + 2x 
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(4.9) 


(4.10) 


(4.11) 


where we consider L = a and N(y(z)) = °D'/2y(x) + (y(z))? in equation (4.11) 


Applying L^! on both sides of (4.11), we get 
y(x) = (y(0) + a?) + L^ [N(y(z))] 


Let us choose, 
volz) = y(0) +27 21 a? 


Then after using (4.6) and (4.7),we get 
Yi+1(x) = L [Ai(yo, yr ,Yi)] > for i = 0, 1, 2, dns 


where A;'s are given by (4.8) or (4.8a). 
Here L-! = J, integral operator and Ag = °D"? yo() + y¿(x) 


yi(z) =I DTS +I EG] 


For = 1 in (4.13): 


yala) = L^! a = dii 1(2) — y1(0) + 21(y0y1)] 
ph ms E 


5/2 i75" 6^ ' 231[7/2 45 


For i — 2, from (4.13) 


13 4 


| 
2 


y3(x) = I! (yo(a)) + 21(yoyz) + L(y?) 


2 4 
eae : z724 + = go? 
2." [772 [9/2 2 21192 


pd uc E. 


(4.12) 


(4.12a) 


(4.13) 


8 


20 
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167 2 1009 139 
5 6 13/2 , 199. 7 
Ter E "a Ze "UC N^ "asi 
1 
98680 gy 17 9 d RS 
17017 x 15|7/2 
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For i = 3, from (4.13) 


mia) = I"? (y3) + 21(yoya) + 21(yiyo) 


a 26 6 z2 13 16 
_ ARA ee o A I 
7z 7t em 9 x R mm^ 


2516 11/2 ER 4 Joan 
6 


+. — 
2315872 ^ 240 ' 5[5/2 [7/2 
For i = 4, from (4.13) 


ys(z) = I"? (y4) + I2yoya + 2y1Y3 + y3) 


z? 8 7 13 1 247 x 16 
A A E UMS A 
3 * g^ + (5 to): ET A 


Similarly, 


(a) if. CA A ( er [5 , 16 
qx) = HH 
is 19/2 Amy 5 8 9x) [mmn 


me fe 79/2 
11/2  9[572 
4 5 56 19 608 
l =1 DO oq? 954.234. 75 45/2 | 7743 7/2 
yADM(2) qe Ge TAT + T Iert 105 z^ 


(ED a- P 1024 ) 9/2 + 
= = T E eee 
12 4m 135./m — 283573/2 


(C) Homotopy Analysis Method (HAM). We consider the same differential 


equation 
Dy(x) — S Data) — (y(z))? = 2x (4.14) 


with initial condition y(0) — 1. 


or, N[y(z)] = Date — *D'?[y(z)] — y? (x) — 2z = 0 (4.15) 
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Using HAM we get, 


Llym(1) — ^mym-1(2)] = REm[Ym-1(1)] (4.16) 
O for m=1 
where Am = and 
1 for m>2 
Rm [Ym—1, m2; 92 Y1, Yo] = z i , 9e NT oi] (4.17) 
(m — 1)! | Ox - 
Here 
N[é(z, q)) = Dé(z; q) - “D'Pg(x, q) — {o(a, a)? — 2x (4.18) 


Expanding in Taylor's Series w.r. to q, 


m=0 q" q=0 
= NT. _ 1 (amd 
= yo(z) + Za ymq™” where yg — — PO (4.19) 
When q = 0, dir 0) = yo(x) and for q = 1, 
p(2,1) =o yi (4.20) 
i=1 
Now, 
T 2 
lola; yy = E vi 
m=0 
= Um (X Um 
m=0 m=0 
co m-1 
= 3. (ym-1-5 yj) à | (4.21) 
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Now, 
gm-1 5 m-1 
gmat UG Doo = A. (m - DlYm-1-5 H 
q j-0 
mo. OU NEUE oc cag (4.22) 
"E ol az m-1-j Yj 
(m — 1)! üqm-1 Ñ 5-0 


From (4.17), from m = 1,2,3,--- 


RmlYo, Y1, ET EN 


m-—1 
"Lg Laeta] 


1 gala 
CA 


q=0 


dë ep1/2 


1 ame 
(m — 1)! V óq- T" 


m-1 jii 
lc aere ur NI E x J MES 
For m 1: 
Run, Vis Y2; *** Men 2. en Al 
= Dym—1 — DP ym-1 — b» UVm-1-j Yj — (1 — Am) (2x) (4.24) 
j=0 
From (4.14): (Ham’s iteration scheme) 
Llyml1) — AmYm-1(2)) 
= h dim —°D Pym — E deu — (1 Ag) 22) (4.25) 
j=0 


We consider linear operator L as the highest order derivative of the differential equation 
(4.14) i.e., L=# and L1 =I. 
Initial condition y(0) — 1 


= yo(0) + y1(0) + ya(0) +---=1 (4.26) 
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Assume yo(x) = 1, then from (4.25), yo(0) = 1 and y;(0) = 0 fori > 1 
From (4.25) for m= 1 


Liya(a)] = h [Dyo —D'Pyo — (yo(2))? - 22] 


Integrating, y, (1) = —h(x + 2?) 
From (4.25), for m > 2 


Ym = Vel) +h 


m-—1 
IDYm-1(2) — IP? ys, tal - I | (Um—j 2] 


j=0 
For m — 2, 
4 16 2 
e Rh +h 2) ua | 3/2 4 g? 5/2, 4 d 
ya (1) (14+h)la+2%) + OC +r to B 
For m = 3, 
: 32 gt 
= (1 A — 973 5/2, 2 7/2, E 
v(e) = (Lb AYE lie (Sue t 
3 4 5 
"TES x z 
" (5 N 2 N a 
For m 4, 
336 3 
= (Te 4 sms ht 5/2 3 7/2 4 
ya(x) = (1+h)[ |+ loe + * 105 yr + 52 
48 oni 5 
taa Trc 
For m = 5, 


z? 8 " 13 1 247 
Im DD t /2 PM c cem dae. mr 09 
Y5 = + (ae)? 3 8 * A(5/2y2 "C gra 


503 8 3032 
+( + E H2 p... 


1155[5/2 ^ 
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For m=6 
z E i bees 
For m = 7, 
yr = e Se mee + 
For h = —1 and a = 1/2 
yo(z) = 1 
yi(z) = +0 
ya(1) = est? +r? + ae = 
y3(x) = E + B wl? y 7 ma aa + Zei + a 
val) = DE Ee? E zy j^ S Geen Hey 
ws) = z+ wei? " N x) K EC SÉ 
yo(x) = NE Ay, qa + 
yr(x) = (s " x) "m 


4 
Unam = lt S ae 


oe) Ic , 1024 ) 9/2 
a bes Gë B ... 
5 135/71 | 283513/2 
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Fig. Graph for linear PDE 
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Fig. Graph for Non-linear PDE 


Conclusion. For comparison of solutions obtained by HAM, ADM and PSM, the 
graphs of the solutions obtained by these methods both for linear and nonlinear PDE. 
It has been found that the graphs corresponding to the solutions upto certain number 
of terms of the above mentioned methods are exactly the same showing that each of 
the different methods has strong base for accuracy of the results, although the question 
of convergency may occur. It can be shown that convergence criterion can be proved 
in Banach space for these methods. 
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